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摘 要 : 种 群 多 样 性 对 微 种 群 教 与 学 优化 算法 的 性 能 有 极 大 影响 。 为 进一步 提高 其 性 能 ， Du s 

样 性 的 微 种 群 教 与 学 优化 算法 (MTLBO-GLD)。 该 算法 从 基因 水 平 上 对 种 群 多样 性 进行 监测 ; 并 使 用 混沌 搜索 和 余 

函数 分 阶段 进行 扰动 以 增加 种 群 多 样 性 。 所 提 算 法 与 8 种 元 启发 式 算法 (4 种 微 种 群 算法 和 4 种 非 微 种 群 算 法 ) 在 
个 测试 函数 上 进行 性 能 比较 。 实 验 结果 表明 ，MTLBO-GLD 的 整体 性 能 要 显著 好 于 其 他 8 种 所 比较 的 算法 。 
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Abstract: The performance of micro-population Teaching-learning-based optimization (TLBO) algorithm is mainly depended 
on the population diversity. In order to further improve its performance, this paper proposed a Micro-population teaching- 
learning-based optimization based on gene level diversity (MTLBO-GLD) . In the proposed algorithm, the population 
diversity is monitored at the genetic level, and the chaos search and the cosine function are used to improve the population 
diversity in different evolutionary stages. The proposed algorithm is compared with eight meta-heuristic algorithms (including 
four micro-population algorithms and four non-micro-population algorithms) on 13 test functions. The experimental results 
show that the overall performance of the MTLBO-GLD is significantly better than that of compared algorithms. 
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0 引言 DESP 算法 可 以 较 快 寻 到 最 优 解 。Salehinejad 等 0 提出 一 种 
E 基于 向 量化 随机 突变 因子 的 差分 进化 算法 (micro-differential 

传统 数学 方法 在 求解 各 类 复杂 优化 问题 时 存在 诸多 不 足 : evolution with vectorized random mutation factor, MDEVM). 
因此 ， 各 类 智能 优化 算法 被 不 断 提 出 。 如 ， 遗 传 算法 (genetic 该 算法 将 突变 因子 向 量 随 机 化 ， 即 给 种 群 中 的 每 个 个 体 随 机 


© algorithm, GA) Ul, SUE SZiX(antcolony optimization, ACO) 分配 一 个 差分 突变 因子 F， 从 而 增加 种 群 的 多 样 性 。 实 验 结 
| 小 、 粒 子 群 算法 (particle swarm optimization，PSO)JB]、 差 分 ; 果 表 明 ， 与 DE 算法 相 比 ，MDEVM 算法 在 搜索 过 程 中 能 够 
— 化 算法 (differential evolution, DE) M4, 较 好 地 增强 种 群 的 多 样 性 , 同时 具有 较 高 的 勘探 能 力 。Brown 


教 与 学 优化 算法 (teaching-learning-based optimization, fij 等 [9 提出 一 种 使 用 小 种 群 的 自 适应 差分 进化 算法 (adaptive 
称 TLBO) 是 Rao 等 人 BSF 2011 年 提出 的 一 种 新 型 智能 优 differential evolution with a small population, uJADE); 该 算法 
化 算法 。 其 原理 主要 是 模拟 学 生 学 习 的 过 程 ， 该 算法 分 为 两 的 改进 是 基于 RerJADE 算法 031。 所 提 算 法 使 用 一 种 新 的 变 
个 阶段 : a) 学 生 在 课堂 上 上 听 老 师 上 课 的 学 习 过 程 ; b) 学 生 课 异 算 子 以 提高 小 种 群 的 全 局 搜索 能 力 ， 并 可 以 加 速算 法 收敛 


T 


的 自我 学 习 过 程 。 由 于 TLBO 算法 实现 简单 、 便 于 理解， 速度 。 此 外 ，hJADE 还 使 用 隔 代 更 新 ucr 和 Ar 的 方式 来 实现 
其 控制 参数 较 少 ， 故 受到 众多 学 者 的 重视 和 研究 ， 并 用 其 解 DERRAT F 和 交叉 因子 CR 的 自 适应 调整 ， 从 而 减少 在 
决 了 许多 实际 问题 。 例如 机 械 优化 设计 器、 热 交换 器 优化 中 、 ”使 用 小 种 群 时 FRI CR 值 迅速 减 小 的 发 生 。 研 究 表明 ， 尽 管 
热电 制冷 器 优化 加、 平面 钢 框架 的 设计 优化 四 等 。 /JADE 算法 使 用 小 种 群 ， 但 从 所 选 测试 函数 来 看 ， 其 表 


一 般 来 说 ，TLBO 算法 及 其 他 智能 优化 算法 都 会 使 用 较 现 比 一 些 使 用 常规 种 群 规模 的 DE 算法 更 好 。Salehinejad 等 
大 的 种 群 规模 来 保持 种 群 多 样 性 以 提高 搜索 性 能 ， 但 在 计算 [提出 一 种 基于 反 疝 学 习 的 集成 差分 进化 算法 (opposition- 
资源 受 限 情况 下 (如 内 存 小 )， 较 大 的 种 群 规模 会 影响 算法 运 based ensemble micro-differential evolution，OEMDE)。 该 算法 
行 。 为 解决 以 上 问题 ， 基 于 微 种 群 的 智能 优化 算法 被 不 断 提 不 仅 使 用 反 向 学 习 生 成 与 当前 个 体 (候选 解 ) 相 反 的 个 体 以 提 
出 。 例 如 ，Ren 等 00 提 出 一 种 小 种 群 的 差分 进化 算法 高 算法 的 全 局 搜索 能 力 ; 在 每 一 代 ， 各 个 个 体 还 随机 选择 一 
(Differential evolution using smaller population，DESP)。 该 算 个 突变 策略 以 提高 种 群 的 多 样 性 。 结 果 表 明 ， 与 其 他 Micro- 
法 不 仅 提 出 一 种 新 的 DE 变异 策略 用 于 提高 全 局 搜索 能 力 ， DE 算法 相 比 ，OEMDE 算法 具有 更 好 的 性 能 。 

还 对 扰动 强度 进行 自 适 应 调整 。 实 验 表 明 , 与 DE 算法 相 比 ， 虽然 众多 学 者 对 微 种 群 优化 算法 进行 了 研究 ， 但 是 他 们 
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的 研究 大 都 基于 个 体 水 平 的 种 群 多 样 性 。 在 实际 种 群 进化 过 
程 中 ， 如 果 某 个 或 某 几 个 诀 策 变量 过 早 收敛 ， 那 么 就 会 对 个 
体 的 进化 造成 极 大 影响 ， 进 而 影响 整个 种 群 进化 。 基 于 文献 
[15] 和 [16] 的 研究 , 本 研究 提出 一 种 基于 基因 水 平 多 样 性 的 微 
种 群 教 与 学 优化 算法 (micro-population teaching-learning- 


在 所 提 算 法 中 ， 首 先 利用 基因 水 平 检测 方法 来 对 微 种 群 算法 
的 种 群 多 样 性 进行 维度 检测 ， 然 后 使 用 混沌 搜索 和 余弦 函数 

行 分 阶段 扰动 以 增强 其 多 样 性 。 最 后 ， 将 所 提 
算法 与 其 他 8 种 算法 在 13 个 测试 函数 上 进行 性 能 比较 。 实 
«8j, MTLBO-GLD 不 仅 具 有 较 快 的 收敛 速度 ， 还 具 
有 较 强 全 局 搜索 能 


1 ”相关 算法 介绍 


1.1 标准 教 与 学 优化 算法 
标准 TLBO 算法 56 主 要 包含 两 个 阶段 , 老师 上 课 的 教学 
过 程 和 学 生 课余 的 自学 过 程 。 具 体 算法 描述 如 下 。 
1.1.1 教学 阶段 
在 教学 阶段 ， 将 适 
阶段 的 教师 ; 
绩 ( 设 为 Xu) 
式 如 下 : 


应 度 值 最 好 的 个 体 ( 设 为 Xn) 作 为 这 个 
教师 希望 通过 教学 过 程 提高 所 有 学 生 的 平均 成 
尽量 靠近 他 自己 的 水 平 。 该 阶段 个 体 的 更 新 公 


Xu m Xiatr(Xo-T, Xy) (1) 
其 中 Xinew 是 经 过 "E 阶段 后 第 i 个 学 生 的 新 位 置 。7F 
round(1t+rand(0,1))， 为 教学 因子 ， 取 1 或 2, r 为 [0,1] 范 围 内 
的 均匀 随机 数 . 若 Xinew 的 适应 度 值 优 于 Xiown， 则 更 新 个 体 ; 
否则 ， 不 更新。 
1.1.2 学 习 阶 段 


based optimization based on gene level diversity, MTLBO-GLD). 
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呈 递 减 趋势 。 因 此 ， 它 不 仅 适 用 于 种 群 的 局 部 扰动 ， 而 且 能 


使 算法 的 进化 后 期 注重 局 部 搜索 。 


N 


根据 文献 [19] 的 研究 , 本 文 


使 用 该 余弦 因子 来 对 个 体 维 度 进 行 扰 动 ， 式 子 如 下 : 


中 , MaxFES 为 设 定 的 最 大 迭代 数 ; N 表示 正 态 分 布 函数 ; 


FES 为 当前 迭代 数 。xjvesio 为 第 g 代 中 多 样 性 最 好 的 维度 ; 


Xjl, 


2 


微 种 群 算法 的 种 群 规模 较 小 ， 


6 为 第 g 代 中 第 姑 个 维度 扰动 后 更 新 的 值 ,jE (1,2,.…, D) 
基于 基因 水 平 多 样 性 的 微 种 群 教 与 学 优化 算法 


对 于 微 种 群 算法 来 说 ， 保 证 种 群 多 样 性 十 分 关键 。 由 于 
导致 种 群 多 样 性 会 和 着 环境 的 


选择 急剧 下 降 ， 这 对 基于 种 群 进化 的 元 启发 式 算法 来 讲 十 分 


不 利 。 
JR, 
it 会 对 整个 个 体 的 进化 带 来 显著 影响 。 因 
[16] 的 研究 , 本 研究 将 对 种 群 进行 基因 水 平 上 的 监测 ， 
使 用 两 种 不 同 的 方法 来 提高 种 群 多 样 性 。 
1 基因 水 平 多 样 性 监测 


HÀ 


2. 


在 学 习 阶 段 中 ， 学 习 者 主要 通过 课 下 与 同学 交流 来 提高 


计算 第 g 代 个 体 变 量 的 第 7 个 维 


同时 ， 已 有 研究 对 于 种 群 多 样 性 的 评判 大 都 基于 个 体 
如 果 种 群 中 的 个 体 出 现 某 一 维 或 某 几 维 收敛 的 情况 ， 

因此 , 基于 文献 [15] 和 
分 别 


基于 文献 [15]， 本 文采 用 方 di 
度 均 值 和 标准 差 公 式 如 下 : 


(5) 


std; = (6) 


一 般 情 况 下 ， 若 第 7 个 维度 的 多 样 性 比较 差 ， 则 其 std;c 


人 也 就 较 小 ,从 理论 上 来 说 , 当 一 个 维度 的 stdic 值 等 于 0 时 ， 


正常 情况 下 ， 元 启发 式 算法 


自己 的 成 绩 。 个 体 的 更 新 公式 如 下 : 度 的 stidic 值 都 是 一 个 逐渐 变 小 的 趋势 。 因此 ,所 提 算 法 先 从 
X as tr (Ex) EF )«rtx) dior i na du 
[agir quid e cedit (2) SIRURSRA UE ESEEREM A AER T ROM THE 
要 。 因此， 本 研究 采取 每 5 代 的 间隔 来 进行 多 样 化 检测 (代数 
其 中 ，p 关 1，p, 1 E(1, 2, .., NP); 有 XY) 为 适应 度 函 数值 。 若 。 ”的 敏感 性 分 析 的 实验 结果 见 3.4.2 节 )。 
Xpnew 的 适应 度 值 优 于 Na 则 更 新 个 体 ， 否 则 ， 不 更 新 。 2.2 多 样 性 的 提高 
1.1.3 标准 TLBO 算法 的 实现 步骤 研究 表明 stds 的 值 降 为 0 需要 很 长 的 时 间 05， 故 本 看 
a) 初始 化 种 群 ， 设 定 维度 D、 种 群 规 模 NP 和 最 大 评价 RH orc 来 识别 种 群 多 样 性 开始 变 小 的 时 刻 ， 式 子 如 下 : 
次 数 。 " -(- D- FES ) l 0D 
b) 教学 阶段 :班级 中 的 每 个 学 员 根 据 式 (1) 进 行 更 新 个 "t MaxFES) | 
体 ， 并 进行 个 体 选 择 。 其 中 
c) 学 习 阶 段 : 通过 式 (2) 对 个 体 进行 更 新 , 并 进行 个 体 选择 。 ajo =max (p06,1—p6) (8) 
d) 判定 程序 是 否 达 到 终止 条 件 , 若 没 有 达到 ， 则 转 至 步 ha We (9) 
Bb); 否则 ， 继 续 执 行 。 upjo - mis 
e) 输出 最 优 解 。 low, = min(m; cX) (10) 
1.2 混沌 映射 up; 7 max(m;o.Xy.;) (11) 
混沌 现象 是 非 线性 动力 系统 的 固有 特性 ， 在 没有 任何 其 中 ，xiowj 和 xy 分别 表示 第 7 维度 预 设 的 上 下 界 。 
种 随机 性 因子 的 作用 下 也 能 够 发 生 类 随机 的 运动 行为 ， 这 种 式 (7) 可 以 知道 orc 的 值 会 随 着 迭代 次 数 的 增加 慢 慢 变 
运动 行为 没有 周期 的 特性 不 具有 收敛 特征 并 且 对 于 初始 化 的 小。 一 般 情 况 下 ， 算 法 种 群 的 多 样 性 也 会 在 迭代 过 程 种 慢 慢 
取 值 有 很 强 的 依赖 性 L071。 混沌 模型 还 有 一 个 重要 的 特征 是 遍 变 差 ， 所 以 本 研究 用 oj6 来 判断 种 群 各 维度 的 多 样 性 。 从 式 
历 性 ， 所 以 用 来 对 个 体 进行 随机 次 数 的 扰动 使 其 跳出 局 部 最 。 (8)~(11) 可 以 看 出 ,qj 主要 为 0c 提供 了 一 些 该 维度 的 信息 ， 


8。 借 鉴 文献 [18] 中 的 混沌 扰动 策略 ， 本 研究 将 利用 
Logistics 搜索 策略 来 对 种 群 进行 基因 水 平 上 的 扰动 , 公式 如 下 : 
Meow = Xa © 40 (1— Xu) (3) 
其 中 ,xnew 是 经 过 混沌 映射 之 后 产生 的 新 变量 ; 4 为 控制 参数 ， 
取 值 范围 在 [0,4] 之 间 ， 且 当 上 =4 Wf, Logistics 映射 是 属于 完 
全 混沌 的 状态 。 

1.3 余弦 扰动 

余弦 函数 


有 比较 强 的 振荡 性 ， 且 在 (0，z/2) 区 间 中 值 


以 便 引 导 ojc 


法 进化 后 期 ， 种 群 的 多 样 性 
会 强调 局 部 搜索 ， 
敛 到 最 优 。 


更 好 的 判断 该 维度 的 多 样 性 。 
此 外 ， 本 文 给 了 一 个 下 限 值 。 一 般 而 言 ， 在 算 
很 小 ， 但 考虑 到 算法 在 进化 后 期 
故 不 再 对 种 群 采取 扰动 措施 ， 直 到 算法 收 
在 本 研究 中 ，ojc 的 最 小 值 设 为 105( 敏 感性 分 析 


ojc 设置 了 


的 实验 结果 见 3.4.3 节 )。 


对 


如 果 监 测 到 某 个 变量 的 多 样 性 丢失 ， 那 么 采用 以 下 方法 


其 进行 扰动 来 提升 多 样 性 。 其 步骤 如 下 : 
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a) 若 10? « std;o ojc, 则 利用 式 (3) 对 该 维度 进行 混沌 扰动 ; ” 平 设 为 5%。 表 中 最 优 结果 加 粗 表 示 。 


b) Xr oa € 10? H.std;c € oi6 , 则 利用 式 (4) 对 该 维度 进 3.1 与 微 种 群 算法 的 比较 
行 余弦 扰动 。 在 这 部 分 中 ， 所 提 算 法 跟 MDEVM, DESP, uJADE X 
2.8 MTLBO-GLD 算法 的 实现 步骤 OBL-uTLBO 4 种 微 种 群 算法 进行 比较 ,所 有 算法 的 种 群 规模 
a) 初始 化 ， 设 定 微 种 群 规模 NP、 维 度 D、 最 大 评价 次 ” 均 设 定 为 8， 最 大 适应 度 函 数 评价 次 数 为 300 7j. 对 于 13 个 
数 ; 初始 化 种 群 ; 测试 函数 的 求解 精度 有 以 下 规定 , 函数 万 的 求解 精度 为 102， 
b) 教学 阶段 :根据 式 (1) 对 个 体 进行 更 新 , 并 对 个 体 进行 选择 ; ER 12 个 测试 函数 的 求解 精度 为 103. 若 算法 在 达到 设 定 的 
c) 学 习 阶 段 : 根据 式 (2) 对 个 体 进行 更 新 , 并 对 个 体 进行 选择 ; 终止 条 件 时 ， 所 求解 还 未 到 达到 规定 的 精度 ， 则 表示 该 算法 
d) 每 隔 5 代 根 据 式 (6) 对 种 群 中 的 每 个 变量 进行 多 样 性 ”在 此 函数 上 求解 失败 ,结果 以 “-” 来 表示 .仿真 结果 如 表 1 所 示 。 
仿 测 操作 ， 并 计算 出 每 个 维度 的 标准 差 ， 然后 根据 式 (7) 计 算 表 1 可 知 ，MTLBO-GLD 算法 对 于 所 有 测试 函数 均 能 
该 维度 相应 的 orc 值 ， 将 两 值 进行 比较 ; 在 终止 条 件 范 围 内 寻 到 最 优 解 ， 且 在 各 个 测试 函数 上 的 求解 
e) 根据 步骤 d) 进行 分 阶段 扰动 : 若 103< stdic< oe, 结果 都 显著 好 于 所 比较 算法 ， 尤 其 在 fi. fo. fs. fa. fe. fio. 


则 利用 式 (3) 对 该 维度 进行 混沌 扰动 ; Æ oe < 103 且 stdic< — 万 这 几 个 测试 函数 上 寻 优 性 能 非常 好 。 其 主要 原因 在 于 ， 对 
ojc , 则 利用 式 (4) 对 该 维度 进行 余弦 扰动 ; 于 MDEVM、DESP 及 HJADE 这 三 个 算法 来 说 ， 虽 然 DE 算 
f) 完成 步骤 di) 和 e 后 ， 对 个 体 进行 选择 ; 法 的 全 局 搜索 能 力 要 好 于 TLBO 算法 , 但 所 提 算 法 通过 从 基 
g) 判定 算法 是 否 达 到 六 mE MAR I 若 结论 是 肯定 ， 寻 水 平 上 对 种 群 进行 监测 ， 并 采用 混沌 与 余弦 函数 分 阶段 来 
则 继续 执行 ， 和 否则 ， 转 回 步 又 b) ; 对 种 群 进行 扰动 。 对 于 OBL-uTLBO 算法 来 说 ， 虽 然 该 算法 
h) 输出 结 使 用 反 向 学 习 来 提高 其 全 局 搜索 能 力 ， 但 并 未 从 基因 水 平 上 

对 种 群 多 样 性 进行 监测 。 因 此 ， 相 较 于 MDEVM、DESP、 


3 (RA HADE 及 OBL-RTLBO 4 种 微 种 群 算 法 ， 所 提 算 法 在 该 函数 
为 测试 所 提 算 法 性 能 ， 本 文选 取 文 献 [20] 中 的 13 个 测试 测试 集 上 的 表现 更 好 。 
函数 ， 并 将 测试 函数 的 维度 设 为 30。 此 外 ， 本 文 还 使 用 非 参数 统计 分 析 方 法 来 对 几 种 算法 所 
将 所 提 算 法 分 别 与 TLBOP 6, DESPU?, MDEVMU!, 求解 的 平均 值 进行 深入 分 析 ， 结 果 见 图 1 和 表 2。 从 图 1 中 
HJADE02]、Rcr-JADE-s403]、FSADEPID、TLBO-CSWLP2 以 及 可知 MTLBO-GLD 算法 的 寻 优 效率 是 最 高 的 , 比 其 中 表现 较 
OBL-AITLBOP3 等 算法 进行 对 比 。 针 对 各 种 算法 的 种 群 规 模 。 好 的 OBL-kTLBO 算法 略 胜 一 筹 ,更 是 显著 高 于 其 他 几 个 比 
会 在 以 下 部 分 进行 说 明 。 所 有 算法 均 在 13 个 测试 函数 上 独 较 算 法 。 另 外 ， 从 表 2 也 可 以 看 出 所 提 算 法 的 寻 优 效率 要 显 
30 次 ， 所 得 结果 取 30 次 的 平均 值 。 同 时 ， 为 保证 所 著 优 于 所 比较 的 几 个 算法 。 
导 结论 的 可 靠 性 ， 仿 真 结果 采用 Friedman""!, Bonferroni- 综 上 所 述 ，MTLBO-GLD 算法 在 基因 水 平 上 对 种 群 多 样 
a Holm 和 HochbergP9g 检 验 进行 统计 分 析 。 显 著 性 水 性 进行 监控 ， 能 够 得 到 更 高 的 求解 精度 。 
del 微 种 群 算法 的 实验 结果 
Tab. 1 Experimental results of Micro-population algorithms 
函数 MDEVM DESP WJADE OBL-pTLBO MTLBO-GLD 


S 


2.8E+04 1.9E+04 2.2E+04 2.6E+02 0.0E--00 
fi 3.3E+03 8.0E+02 7.8E+02 8.9E+00 0.0E--00 
3.1E«04 i 3.7E+04 3.9E+02 0.0E--00 

f 4.7E+03 s 12E403 1.1E+01 0.0E--00 
1.4E+05 E 1.6E+05 4.3E+02 0.0E--00 

n 1.5E+04 * 7.7E+03 5.3E+01 0.0E+00 
2.2E+05 4.3E+02 0.0E--00 

f : " 1.6E+04 1.2E+01 0.0E+00 
2.3E+06 2.1E+05 1.8E+05 2.SE+01 

f 4.0E+05 - 5.5E+04 4.3E+04 1.2E+00 
1.8E+04 à 1.2E+04 1.2E+02 0.0E--00 

f 6.3E+03 : 3.8E+03 1.0E+01 0.0E+00 
8.9E+05 z 2.3E+05 2.1E+03 2.5E-06 

f 标准 方 7.1E+05 - 1.6E+05 8.1E+02 1.8E-06 
平均 1.0E+05 1.9E+05 3.8E+03 

f - : 5.1E+03 6.7E+03 1.4E+03 
. i 1.2E+05 2.7E+02 0.0E+00 

f - B 6.8E+03 1.3E+01 0.0E+00 
" 3.8E404 3.9E402 0.0E--00 

fo - : 5.0E+03 1.0E+01 0.0E+00 
2.6E404 3.1E404 4.6E+04 2.8E+02 0.0E--00 

s 3.1E+03 1.7E+03 3.7E+04 1.3E+01 0.0E--00 
6.6E+04 1.8E+06 3.8E+04 4.1E+03 6.9E-05 

fs 6.5E404 0.0E+00 7.5E+03 7.8E+02 1.9E-04 
7.0E+04 i 3.2E+04 2.5E+05 1.1E-02 

fe 6.4E404 : 1.8E+04 4.8E+05 1.2E-02 
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5 4.5769 Hochberg’s procedures for Micro-population algorithms 
n 3.9615 | MTLBO-GLD vs. z  Unadjusted p Bonferroni-Dunn p — Holm p Hochberg p 
aa 3.1923 | DESP 5.77 8.04E-09 3.22E-08 3.22E-08 3.22E-08 
È 2.2692 MDEVM 4.78 1.79E-06 7.18E-06 5.38E-06 5.38E-06 
0 | uJADE 3.53 4.08E-04 1.63E-03 8.16E-04 8.16E-04 
ir M 1 OBL-HTLBO 2.05 4.07E-02 1.63E-01 4.07E-02 4.07E-02 
° MTLBO-GLD OBL-hTLBO ur MDEVM DESP a 与 非 微 种 群 算法 的 比较 
在 这 部 分 中 ， 所 提 算 法 跟 FSADE、Rcr-JADE-s4、TLBO 
图 1 Friedman 测试 在 —— fI TLBO-CSWL ix 4 种 非 微 种 群 算法 的 性 能 进行 比较 。 期 户 
Fig. 1 Ranking obtained by Friedman's test for micro-population algorithms 求解 精度 同 3.1 部 分 。 四 个 算法 的 种 群 规模 设 定 
表 2 Bonferroni-Dunn、Holm 和 Hochberg 在 微 种 群 算法 上 所 得 p fü 依次 为 50、100、30 和 40， 算 法 参数 设置 均 与 文献 [12] 
Tab.2 P-values obtained by Bonferroni-Dunn's, Holm's, and 一 致 。 结 果 见 表 3 。 


表 3 所 提 算 法 与 非 微 种群 算 法 的 结果 


Tab.3 Experimental results of the proposed algorithm and non-Micro-population algorithms 


FSADE Ro-JADE-s4 TLBO TLBO-CSWL MTLBO-GLD 
函数 (NP=50) (NP=100) (NP=30) (NP=40) (NP=8) 
均值 3.5E+04 2.1E+04 4.3E+03 1.9E+03 0.0E+00 
fi 标准 方差 5.3E+02 4.9E+02 9.1E+02 4.2E+01 0.0E+00 
平均 值 4.5E404 3.3E+04 7.6E+03 2.8E+03 0.0E+00 
f 标准 方差 5.7E*02 7.6E*02 1.3E*02 4.6E+01 0.0E--00 
平均 值 2.6E+06 6.7E+04 1.7E+04 2.9E+03 0.0E+00 
f 标准 方差 4.8E+05 3.6E+03 8.6E+02 1.6E+02 0.0E+00 
平均 值 2.2E+05 1.0E+05 7.2E+03 3.2E+03 0.0E+00 
^ 标准 方差 4.2E+03 3.5E*04 1.5E*02 5.1E*01 0.0E--00 
平均 值 2.8E+06 7.8E+04 - 4.1E+05 2.SE+01 
f uis 
标准 方差 6.1E+04 2.7E+03 - 5.8E+04 1.2E+00 
平均 值 1.3E+04 8.8E+03 1.7E+03 1.1E+03 0.0E+00 
f 标准 方差 4.7E+02 3.6E+02 2.2E+02 2.1E+02 0.0E+00 
平均 值 1.1E+05 2.2E+04 - 2.2E+04 2.SE-06 
fi sS 
标准 方差 2.6E+04 9.8E+03 - 6.4E+03 1.8E-06 
平均 值 4.4E+04 7.0E+04 - 4.5E*05 3.8E+03 
f re 
标准 方差 1.8E+03 2.8E+03 z 8.2E+03 1.4E+03 
平均 值 8.9E+04 9.1E+04 7.9E+03 0.0E+00 
^ 标准 方差 4.8E+03 1.3E+03 8.6E+03 0.0E--00 
平均 值 4.8E+04 3.1E+04 2.3E+04 2.8E+03 0.0E+00 
n 标准 方差 5.0E+02 6.4E+02 1.3E+02 4.4E+01 0.0E+00 
平均 值 3.9E+04 2.2E+04 4.1E+03 2.1E+03 0.0E+00 
fi 标准 方差 3.7E+03 6.4E+02 3.9E+02 5.8E+02 0.0E+00 
平均 值 4.0E+04 1.9E+04 - 5.0E+04 6.9E-05 
fia mn 
标准 方差 1.0E+03 5.5E+02 - 2.5E+03 1.9E-04 
平均 值 4.1E+04 2.1E+04 - 3.8E+05 1.1E-02 
fis 1.1E+03 4.9E+02 6.1E+04 1.2E-02 
从 表 3 中 可 以 看 出 ，MTLBO-GLD 算法 在 13 个 测试 函 羽 水 平 上 对 种 群 进行 监控 是 一 种 十 分 有 效 的 方法 。 
数 上 所 得 结果 均 好 于 其 他 几 个 非 微 种 群 算法 。 因 此 ， 所 提 算 5 r 
法 虽然 使 用 的 是 较 小 规模 的 种 群 ， 但 它 仍 然 表现 出 非常 好 的 外 39231 i) | 
寻 优 性 能 。 原 因 主要 有 以 下 两 点 ， 首 先 ， 所 提 算 法 从 基因 水 | um | 
平 上 对 种 群 的 多 样 性 进行 监测 ， 从 而 可 以 从 种 群 中 提取 更 多 E rs 
效 信息 ; 其 次 ， 采 取 了 分 阶段 对 种 群 进行 扰动 ， 从 而 可 以 TM 
保持 算法 种 群 的 多 样 性 。 因 此 ， 算 法 的 寻 优 性 能 得 到 了 很 大 E E 
的 提升 。 9 MTLBO-GLD TLBO-CSWL nu s4 TLBO FSADE 
为 使 实验 结论 更 加 可 靠 ， 还 使 用 非 参数 的 统计 方法 来 对 
所 有 算法 的 均值 进行 分 析 。 统 计 分 析 结果 见 图 2 和 表 4。 由 图 2 Friedman 测试 在 MTLBO-GLD m — 


图 2 可 知 ， 与 使 用 常规 规模 种 群 的 算法 相 比 ，MTLBO-GLD Fig.2 Ranking obtained by Friedman's test for MTLBO-GLD and 
的 寻 优 性 能 排名 第 一 。 从 表 4 可 以 看 出 ，MTLBO-GLD 算法 non-micro-population algorithms 
的 求解 性 能 要 显著 性 好 于 所 比较 的 所 有 非 微 种 群 算法 。 3.3 MTLBO-GLD 5 OBL- u TLBO, TLBO-CSWL 的 比较 


综 上 所 述 ， 所 提 算 法 虽然 为 微 种 群 算法 ， 但 其 寻 优 性 能 为 了 进一步 验证 在 资源 受 限 情况 下 的 MTLBO-GLD 算 
不 比 常规 种 群 规模 的 算法 差 ， 反 而 更 优 。 这 恰好 证 明 从 基 ”法 性 能 ， 这 里 本 文 将 对 所 提 算 法 与 前 面 两 部 分 中 表现 较 好 的 
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OBL- u TLBO 和 TLBO-CSWL 进行 性 能 比较 。 
的 是 ,在 本 实验 中 ,最 大 的 适应 度 函 数 评 价 次 数 设 为 30,000。 
时 ，OBL-uTLBO fll TLBO-CSWL 的 种 群 规模 分 别 设 为 8 
和 40。 所 得 结果 见 表 5. NUK 5 可 以 看 出 MTLBO-GLD 算法 
TE fi. fo» fs fas fs. fe. fios fui. fis ESSSIBSfEEE TLBO-CSWL 
算法 的 解 要 好 很 多 。 而 TLBO-CSWL 只 在 这 个 函数 上 要 显 


与 前 面部 分 不 


I 


Hl 


著 好 于 所 提 算 法 。 结 果 足 以 表明 所 提 算 法 的 性 能 是 好 于 
TLBO-CSWL 算法 。 对 于 OBL-nTLBO， 两 个 算法 在 这 些 测 
试 函数 上 得 到 的 平均 结果 总 体 上 来 说 差不多 ， 但 所 提 算 法 在 
每 个 测试 函数 上 表现 的 比 OBL-nTLBO NE EX Lib, 
MTLBO-GLD 不 仅 具 有 较 快 的 收敛 速度 ， 还 能 够 降低 计算 成 
本 ， 同 时 也 展现 了 较 强 的 全 局 搜索 能 

本 部 分 还 使 用 非 参 数 统计 分 析 方 法 来 对 三 种 算法 的 平均 
值 进行 分 析 ， 其 结果 见 表 6。 由 表 6 可知 MTLBO-GLD 算法 


的 整体 性 能 优 于 TLBO-CSWL 算法 。 另 外 ， 昌 然 MTLBO- 
GLD, TLBO-CSWL, OBL-uTLBO 的 性 能 在 统计 学 意义 上 不 
存在 明显 的 差异 ， 但 是 从 表 6 中 的 结果 可 知 , MTLBO-GLD 
的 整体 性 能 仍 稍微 优 于 OBL-&TLBO 和 TLBO-CSWL 的 。 
K4 Bonferroni-Dunn, Holm 和 Hochberg 

在 MTLBO-GLD 和 非 微 种 群 算法 上 所 得 bp 值 

Tab.4 P-values obtained by Bonferroni-Dunn's, Holm's, and Hochberg’s 


procedures for MTLBO-GLD and non-micro-population algorithms 


MTLBO-GLD vs. z Unadjusted p Bonferroni-Dunnp Holmp Hochberg p 


表 6 所 有 比较 算法 的 Bonferroni-Dunn、Holm 以 及 Hochberg 检验 的 p 值 
Tab.6 P-Values obtained by Bonferroni-Dunn's, Holm's, and Hochberg's 


procedures on experimental results with all comparison algorithms 


MTLBO-GLD vs. z  Unadjustedp Bonferroni-Dunnp Holmp Hochberg p 


TLBO-CSWL 1.86 6.24E-02 1.25E-01 
OBL-uTLBO 0.49 6.24E-01 1.00E+00 


1.25E-01 1.25E-01 
6.24E-01 6.24E-01 


3.4 实验 分 析 
3.4.1 混沌 和 余弦 两 种 扰动 的 结果 分 析 

一 般 而 言 ， 算 法 在 前 期 应 偏重 于 全 局 搜索 ， 在 后 期 应 偏 
重 于 局 部 搜索 RP。 在 本 研究 中 ， 所 提 算 法 采用 混沌 搜索 和 余 
弦 扰 动 两 种 方法 进行 分 阶段 增强 种 群 多 样 性 。 为 验证 混沌 


索 和 余弦 扰动 两 种 方法 的 特点 ， 本 部 分 给 出 两 种 方法 在 二 维 
空间 中 的 个 体位 置 变 化 图 。 这 里 以 8 个 粒子 为 例 ， 粒 子 迭 代 
60 次 的 位 置 变化 图 如 图 3 和 4， 从 图 3 和 4 中 可 以 看 出 ， 混 
沌 搜索 的 全 局 搜索 能 力 要 好 于 余弦 函数 。 
100 T T T L7. T 
dH n . o : e 
di 2 : XS . : 
40 * tet T. 9 ~ v. 9e 
setti iba de uus 
tu UN [SU Y PES iuh 
op * ?$* . e ^% 
Lo M ^t we - 4 : 
MEG CEDE 
< 四 e " I E L is . geo . i ° 
60 + r] * e. ë 
80 ? [P hd 


1 1 1 1 1 
0 20 40 60 80 


FSADE 5.0854 3.67E-07 1.47E-06 1.47E-06 1.47E-06 

TLBO 4.7133 2.44E-06 9.75E-06 7.31E-06 7.31E-06 

Ra-JADE-s4 3.7831 1.55E-04 6.20E-04 3.10E-04 3.10E-04 

TLBO-CSWL 2.0225 4.97E-02 0.2182 4.97E-02 4.97E-02 

TLBO-CSWL 2.5427 1.10E-02 4.40E-02 1.10E-02 1.10E-02 
表 5 所 有 比较 算法 在 30 维 函 数 上 所 得 结果 


Tab.5 Experimental results of all compared algorithms on 30D functions 


TLBO-CSWL OBL-uTLBO MTLBO-GLD 
(NP-40) (NP-8) (NP-8) 
F 4.0E-223 0.0E--00 0.0E--00 
f auae 0.0E+00 0.0E--00 0.0E--00 
Hg 9.8E-141 0.0E--00 0.0E--00 
Poo augsx  6.5E-140 0.0E--00 0.0E--00 
平均 值 2.8E-221 0.0E--00 0.0E--00 
PO agg 0.0E+00 0.0E--00 0.0E--00 
3548 5.5E-224 0.0E--00 0.0E--00 
^00 aug 0.0E+00 0.0E+00 0.0E--00 
FH 2.8E+01 2.8E+01 2.6E+01 
f 0 aug 4.1E-01 2.78E-01 3.5E-01 
Pt 9.2E-16 0.0E+00 0.0E+00 
^ 0 BEA 1.5E-15 0.0E+00 0.0E+00 
P 0.0E+00 4.4E-05 2.5E-04 
P — ED%  0.0E+00 3.2E-05 1.4E-04 
平均 值 0.0E+00 6.4E+03 4.0E+03 
^ EI%  0.0E+00 1.3E+03 1.5E+03 
平均 值 0.0E--00 0.0E--00 9.7E-01 
^ —— aux 0.0E+00 0.0E--00 2.5E«00 
FH 4.9E+03 2.1E-16 0.0E+00 
fo aue 1.1E+03 8.5E-16 0.0E--00 
P 0.0E+00 0.0E+00 0.0E+00 
fOO aue 0.0E--00 0.0E--00 0.0E--00 
P 3.2E-219 8.3E-02 1.9E-04 
fo age 0.0E--00 9.7E-02 2.8E-04 
平均 什 2.9E401 2.6E«00 3.1E-03 
f^ — p — 630E-02 5.3E-01 6.2E-03 


图 3 混沌 搜索 
Fig. 3 Chaotic search 
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图 4 余弦 扰动 
Fig. 4 Cosine perturbation 

3.4.2 间隔 代数 值 敏感 性 分 析 

本 部 分 对 间隔 的 代数 值 进行 敏感 性 分 析 。 本 实验 中 ， 种 
群 规模 设置 为 8， 最 大 的 适应 度 函 数 评价 次 数 设 为 30,000, 
MTLBO-GLD 算法 独立 运行 30 次 ;并 选取 13 个 测试 函数 。 
间隔 代数 分 别 设 定 为 0O、5、10、15 和 20， 并 使 用 非 参 数 的 
统计 方法 对 平均 值 进行 排序 分 析 。 各 个 参数 的 排序 结果 如 医 
5 所 示 。 由 图 5 可 知 ， 当 间隔 代数 值 取 S 时 ， 所 提 算 法 的 整 
体 性 能 表现 最 佳 ， 所 以 间隔 的 代数 值 设 为 5。 
3.4.3 oj,G 值 敏感 性 分 析 

该 部 分 对 oio 的 下 限 值 进行 敏感 性 分 析 ， 实 验 设置 与 
3.4.2 相同 。ojc 分 别 设 定 为 1077. 102, 103, 1047 105, 不 
同 参 数 下 的 排序 结果 见 图 6。 由 图 6 可 知 , 在 orc 的 下 限 设 为 
105 时 ， 所 提 算 法 的 整体 性 能 最 佳 。 因 此 ， 在 所 提 算 法 中 ， 
oc 的 下 限 值 设 为 105。 
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图 5 Friedman 测试 在 不 同 代数 间隔 下 所 得 排序 结 


Fig. 5 prd: obtained P Friedman's test for different iteration intervals 
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图 6 Friedman 测试 在 oj,G 取 不 同 下 限 值 下 所 得 排序 结果 
Fig.6 Ranking obtained by Friedman's test for oj,G under different lower limits 


4 ”结束 语 


为 平衡 微 种 群 教 与 学 优化 算法 的 全 局 和 局 部 搜索 能 力 ， 
本 文 提出 一 种 基于 基因 水 平 多 样 性 的 微 种 群 教 与 学 
(MTLBO-GLD) 算 法 。 在 该 算法 中 , 首先 从 基因 水 平 上 对 种 群 
进行 监测 ， 以 防止 某 一 维 或 某 几 维 过 早 收敛 而 对 整个 种 群 搜 
索 带 来 不 利 影响 ， 然 后 使 用 混沌 搜索 和 余弦 函数 分 阶段 实现 
i 动 来 实现 种 群 多 样 性 。 所 提 算 法 跟 JADE, DESP, MDEVM 
OBL-uTLBO 等 微 种 群 算法 和 FSADE、Rcr-JADE-s4、TLBO、 
TLBO-CSWL 等 非 微 种 群 算法 进行 性 能 比较 。 实 验 结果 表明 
MTLBO-GLD 算法 的 整体 性 能 要 优 于 所 有 比较 算法 。 
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